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We measure the quantum speed of the state evolution of the field in a weakly-driven optical cavity QED 
system. To this end, the mode of the electromagnetic field is considered as a quantum system of interest 
with a preferential coupling to a tunable environment: the atoms. By controlling the environment, i.e., 
changing the number of atoms coupled to the optical cavity mode, an environment assisted speed-up is 
realized: the quantum speed of the state re-population in the optical cavity increases with the coupling 
strength between the optical cavity mode and this non-Markovian environment (the number of atoms). 
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Identifying time-optimized processes is a ubiquitous 
goal in virtually all areas of quantum physics, such as quan¬ 
tum communication [1], quantum computation [2], quan¬ 
tum thermodynamics [3], quantum control and feedback 
[4], and quantum metrology [5]. To this end, the notion 
of a quantum speed limit (QSL) has proven to be useful 
and important. The QSL determines the theoretical upper 
bound on the speed of evolution of a quantum system. It 
can be understood as a generalization of the Heisenberg un¬ 
certainty relation for energy and time. It has been derived 
for isolated, uncontrolled systems [6-8], time-dependent 
Hamiltonians [9-13], and more recently for more general 
open system dynamics [14-21]. Although fundamental in 
nature, practical consequences or even experimental appli¬ 
cations of the QSL are still lacking. Nevertheless, achiev¬ 
ing the maximal quantum speed is of high practical rele¬ 
vance, especially in the development of quantum informa¬ 
tion processing devices. 

On the theoretical side, a recent study [14] hinted at 
the possibility of observing speed-ups of the quantum evo¬ 
lution if an open quantum system is subject to environ¬ 
mental changes. Ref. [14] analyzes the dynamics of the 
damped Jaynes-Cummings model, which describes many 
cavity QED systems. These systems in both the interme¬ 
diate and strong coupling regime can exhibit environment- 
assisted evolution [22] - such as non-exponential decay. 

This letter reports an experimental realization of the the¬ 
oretically proposed environment assisted speed-up [14]. To 
this end, we look at the system in an unusual way - as just 
consisting of the cavity field. This allows us to treat the 
atomic number that generates the atomic polarization (the 
off diagonal elements of the atomic master equation) as a 
tunable environment with a range of coupling constants. 
We demonstrate that increasing the interaction of the opti¬ 
cal cavity field with the environment by tuning the number 
of atoms, modifies the time dependent non-classical inten¬ 
sity correlation function, enhancing - speeds-up- the rate 
of evolution of the cavity field in a range with no clear os¬ 


cillations present. 

Our cavity QED system operates in the intermediate cou¬ 
pling regime, where the cavity-atom parameters are of the 
same order: (g, ft, 7 ) / 2 tv = (3.2,4.5, 6.0) MHz. Here g 
denotes the electric dipole interaction strength of an atom 
maximally coupled to the cavity mode (also known as the 
single photon Rabi frequency), and ft and 7 are the de¬ 
cay rates of the cavity and the atom, respectively. For 
weak driving and N atoms in the cavity [23], the effec¬ 
tive dipole coupling between cavity and atomic environ¬ 
ment scales as gy/N « Q VR the vacuum Rabi Oscilla¬ 
tion. The interaction of cavity and atoms is stronger for 
larger N. In contrast to conventional studies of cavity QED 
we make full use of this observation and its tunability by 
using a slow atomic beam: The atomic beam, i.e., a col¬ 
lection of N two-level atoms randomly positioned in the 
cavity mode, with a range of position dependent coupling 
constants g(r, 9 , z), quantified by N eff atoms [24], can be 
understood as a controllable, tailored environment for the 
cavity mode, which we consider as the quantum system (of 
interest). Our atomic beam can produce a wide range of ef¬ 
fective number of atoms in the cavity (W eff = 0.1 -A 30), 
which allows us to investigate the speed of state evolution 
as we change the number of atoms in a controllable way. 

Intuitively, one expects to observe a faster state evolution 
when more atoms are in the cavity; the more subtle effect 
is that there is a non-linear dependence on the number of 
atoms. Conditional measurements of the photons leaving 
the cavity ( g^ (r), second order intensity correlation func¬ 
tion) are ideal to study these environmental effects. In our 
experiment we measure g^ (r) as we change the environ¬ 
ment, the number of atoms. More concretely, we investi¬ 
gate the initial antibunching dynamics as the state returns to 
a steady state before vacuum Rabi oscillations are present 
[24-26]. 

Our measurements clearly show the theoretically pre¬ 
dicted environment-assisted speed-up with its non-linear 
dependence on the number of atoms. This study not only 
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verifies fundamental predictions, but also opens new av¬ 
enues to control the quantum-speed-limited dynamics in 
any cavity QED system. 

Theoretical predictions. We start by discussing a sim¬ 
ple, phenomenological model for our system. To this end, 
we consider the optical cavity to be weakly driven by an ex¬ 
ternal field of strength proportional to e, which generates a 
field coupled to N two-level atoms through the Coopera- 
tivity C — NCi, C\ — g 2 /^. This situation is described 
by the master equation [23] 

pit) =e[a ] - a , p(t)] + g [a f J_ - aJ+,p(t)\ 

+ k (2a p(t) a * — a)a p(t ) — p(t) a^a) 

+ 7/2 ( 2<7 - P(t) °+ “ a + a - P(t) ~ P(t) a + a -) • 

3 = 1 

( 1 ) 

Here a) and a describe the creation and annihilation opera¬ 
tors for the cavity mode, g\ and a J _ correspond to the Pauli 
spin matrices for the jth atom, while J± = J2j=i a ± are 
collective excitations. In the limit of weak driving, e/k, <C 
1, a perturbative treatment is possible. Then, a solution of 
Eq. (1) can be approximated by p(t) ~ || with 

\m) * 100) + A 1 (t) 110) + A 2 (t) |01) + 0{e 2 /k 2 ) 

( 2 ) 

where | nm) state with n photons in cavity and m excited 
atoms. The dynamics of the amplitudes Ai(t) and A 2 (t) 
are described by [23] 

A^t) =-nA^t) + gy/N A 2 (t) + e 
A 2 (t) = --y/2A 2 (t) - gVN A ± (t). 

Adiabatically eliminating the cavity k <7 7 , gy/N one gets 
the Purcell enhancement factor that changes the exponen¬ 
tial decay rate of the atoms from 7 to 7(1 + 2 C) which 
depends linearly oniV. A similar result for the cavity de¬ 
cay rate is found by adiabatically eliminating the atoms 
7 <C n,gy/N where now the decay rate of the cavity 
changes from n to n(l + 2C) again linearly dependent 
on N. Here we want to stress a different regime, where we 
keep both entities in the evolution and so we immediately 
observe that the cavity field couples to N two-level atoms 
through a dipole transition with rate gy/N - an effect not 
indicated by oscillations. Thus, by tuning the number of 
atoms N the interaction strength of the cavity mode and 
the atoms can be varied. In the following, we focus exclu¬ 
sively on the escaping field of the system as we vary the 
number of atoms. Thus, we can monitor the changes in its 
dynamic behavior as we modify this “variable reservoir” - 
showing the potential of engineering dynamical behavior 
by manipulating the reservoir. 

Reference [23] also derived the solution of Eq. (3) and a 
closed expression for the second order correlation function 



FIG. 1. (color online) Calculated correlation function g^ 2 \r) 
(4) accounting for experimental effects for (g,K, 7 )/ 27 r = 
(3.2, 4.5, 6.0) MHz with Q VR /2n = 1.1 MHz (black, continu¬ 
ous line), 2.8 MHz (red dotted line), and 5.2 (blue, dashed line). 
The inset corresponds to a model with maximally coupled atoms 
and optimal experimental conditions. 


g( 2 ) (t) applicable for N atoms, which reads 

9 m (T ) = { 1 + ^ eX p(_r («+!)) 
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where the parameters are defined by, A a/a = 
-2C((2C/(l + 2C- 2 CO), and C( = CV(1 + 7 /2k). 
Furthermore, the vacuum Rabi frequency is, f2y R = 
V((«-7/2)/2) 2 -s 2 iV. 

This master equation model (1) assumes static, maxi¬ 
mally coupled atoms to the cavity field with zero cavity 
and atomic detuning. A refined model simulates fluctua¬ 
tions of the atomic beam number with a Poissonian weight 
and randomly generates the positions of the atoms in the 
cavity mode (which is radially Gaussian with a longitudinal 
standing wave) [27, 28]. This leads to variable dipole cou¬ 
pling strengths and effectively decreases the overall cou¬ 
pling. The atoms present then a range of couplings and 
behave as a reservoir depending on where they are located 
on the mode. The calculation also includes experimental 
effects from cavity and atomic detuning of g 2 {r) is plotted 
in Fig. 1 for the corresponding £Ivr [26, 29]. We observe 
that while the Rabi oscillations are strongly suppressed the 
rate of refilling characterized by the slope at the inflection 
point, becomes a characteristic for the state dynamics re¬ 
turning to the steady state. The atoms, with their different 
couplings constants to the cavity mode behave as an inho- 
mogeneously broadened reservoir, with a variety of cou¬ 
pling constants that nevertheless preserve the non-classical 
anti-bunched nature of the field. All the oscillations start 
anti-bunched but their oscillations are different and averag¬ 
ing them produces dephasing that results in a simple decay 
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rate that will preserve the non-linear y/N dependence. The 
inhomogenous reservoir is crucial for this discussion. The 
inset in Fig 1 shows the result of N maximally coupled 
atoms and optimal experimental conditions. 

Before we continue the discussion with the experimen¬ 
tal results, note that the speed-up proposed here is phys¬ 
ically equivalent to the behavior of the QSL predicted in 
Ref. [14]. However, here we are interested in the dynamics 
of the field, whereas Ref. [14] focused on the atoms. 

Experimental set-up and results. Figure 2 shows the 
general layout of the physical system, which is our new- 
generation optical cavity QED apparatus based on the pre¬ 
vious one described in Ref. [30] operating on the D 2 line 
of Rb. The new system has n sat = y 2 /3g 2 = 1.2 and 
Ci = g 2 /n 7 = 0.38. An imbalanced magneto optical 
trap (MOT) produces a slow atomic beam of 85 Rb atoms 
that travels down and couples to the TEM 0 o mode of our 
0.8 mm long optical cavity. The magnetic field points along 
the —z axis and defines our quantization axis. The MOT 
coils and an arrangement of four permanent magnets gen¬ 
erate this field of about B — — 7.2 G z. 
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FIG. 2. (color online) Simplified diagram of the apparatus for 
measuring autocorrelations with an atomic beam in the presence 
of a magnetic field B. HWP: half wave plate; PBS: polarizing 
beam splitter; BS: beam splitter; APD: avalanche photo-diode. 


The downward pointing laser beam from the MOT pro¬ 
vides optical pumping into the \F = 3, m F = 3) state. We 
verify this by measuring the absorption as a function of fre¬ 
quency, keeping the cavity resonant to the laser frequency. 
We drive the system with H -polarized light (perpendicu¬ 
lar to the B magnetic field). In the frame of the atoms, it 
appears as a combination of two circularly polarized fields, 
while only one is capable of driving with the stretched tran¬ 
sition. 

Experimentally, we increase the number of atoms by 
raising the current of the Rb dispenser, which increases the 
number of atoms in the MOT. A locked 820 nm laser serves 
to keep the cavity on resonance with the atoms and the ex¬ 
citation laser at 780 nm. 

Polarization elements and mode-matching optics prepare 
the driving laser before it enters the cavity. A lens at the 
output collimates the beam and a half-wave plate (HWP) 
aligns the polarization to a calcite PBS which separates the 
H (driven) and V (undriven) modes. The H mode passes 
through a second beam splitter which divides the light be¬ 
tween two avalanche photodiodes (APD) for single or co¬ 


incidence measurements. A series of frequency and spa¬ 
tial filters remove background light. The APD electronic 
output pulses then go to a time-stamp unit in a computer 
which, after post-detection signal processing, produces the 
appropriate correlations. The total detection efficiency of 
the system is 30%. We block part of the laser beam that 
could co-propagate with the atomic beam either with a dark 
spot or with an axicon system that produces a dark region at 
the cavity [31]. The geometry and the detection that we use 
place limitations on the state purity and signal background. 
Nevertheless, we can identify the non-classical feature of 
antibunching and measure the time the system needs to re¬ 
lax back to its steady state. 

We measure QyR^eff) in the cavity QED system by 
transmission spectroscopy for a given the atomic flux. We 
scan the drive frequency while keeping the atoms and the 
cavity resonant [32, 33]. The separation between the two 
peaks is twice Q V r- We observe a clear splitting on the 
transmission spectroscopy while on the correlation func¬ 
tions the corresponding time oscillations are not present 
[26] for the range of parameters that we explore in this 
paper; however we do observe the oscillatory regime for 
higher dispenser currents in the MOT. This is related to the 
different kinds of averaging that take place in the transmis¬ 
sion spectra and in the correlation functions. 
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FIG. 3. (color online) Measured g^ 2 \r) (red circles with only 
statistical error bars) showing antibunching for ^vRV^Aeff = 
8.5 =b 0.1 MHz with fit to an inverted Lorenzian (blue, solid line) 
with the amplitude ao and the HWHM indicated, as well as the 
prediction from the refined model (green, dashed line). 

Figure 3 shows one of our experimental traces of g^ (r). 
We extract the quantum speed from the rate of change of 
the anti-bunching towards the steady state. g 2 (r) reaches 
unity only after a long time, which is a remnant of the tran¬ 
sit time of the atoms through the cavity mode [26]. We 
fit an inverted Loretzian (continuous blue line), with three 
parameters: offset, amplitude a 0 , and the half width, half 
maximum (HWHM). This gives a quantitative measure¬ 
ment of the rate of refilling though the slope by taking the 
ratio a 0 /HWHM, The dashed green line shows the results 
of the refined model (as for Fig. 1) that includes experimen- 
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FIG. 4. (color online) Speed at the inflection point, ao/HWHM, 
(blue circles) as function of f^vR from measured antibunching 
curves. Inset ao/HWHM from the refined model (red dots) as in 
Figs. 1 and 3. 

tal fluctuations on the cavity detunings (±2.5ft), the pres¬ 
ence of a Zeeman sub level 5 MHz away, and the atomic 
spatial distribution. The amplitude is scaled to take into 
account experimental backgrounds. Spurious electronic 
noise makes the fluctuations larger than purely statistical. 

Figure 4 shows the growth of the rate of refilling (blue 
circles) as a function of Q vR , for which we expect a linear 
dependence. The dashed straight line is a fit with slope 
0.24 dt= 0.03 /xs -1 /MHz, with reduced x 2 — 1-67. The 
inset (red dots) shows the result of a simulation with slope 
0.29 =b 0.02 i± s -1 /MHz, of the refined model. 

Our results clearly indicate a trend of enhanced evolution 
speed as a function Q V r « gy/^es- This is in agreement 
with theoretical predictions [14] that observe a speed-up 
as the coupling increases in a cavity QED system and in 
contrast to the dependence of N in the adiabatically elim¬ 
inated cases stated above, that do not show in the non- 
classical correlation function. Additional measurements 
show a similar linear behavior of the slope at the inflection 
point when we have a bunching peak. The speed increases 
linearly as a function of 12 VR . 

Non-Markovian speed-up . In the example studied in 
Ref. [14] the quantum speed-up in an open quantum sys¬ 
tem was attributed to the non-Markovian nature of the en¬ 
vironment. The obvious question remains whether the here 
reported speed-up is also a non-Markovian effect. 

Generally, quantifying non-Markovianity is complicated 
[22, 34, 35]. However, for simple quantum systems such as 
the Jaynes-Cummings model the situation drastically sim¬ 
plifies and a measure of non-Markovianity (J\f) can be in¬ 
troduced [36]. It is defined as 

J\T = max / dt a(t, pi 2 (0)) (5) 

pi,2(0) J a> 0 

where pi j2 (0) are two initial states and 

cr(*, ^ 1 . 2 ( 0 )) = ~trflMt) ~ P2(t)\} . (6) 


It has been seen [36] that for Markovian dynamics all ini¬ 
tial states monotonically converge towards a unique sta¬ 
tionary state. Thus, cr(£, pi, 2 (0)) (6) is strictly negative 
and J\f = 0 (5). Non-Markovian dynamics are charac¬ 
terized by an information backflow from the environment, 
and the convergence of p(t) towards the stationary state is 
accompanied by oscillations. Hence, cr(£, pi j2 (0)) (6) can 
become positive, which amounts to finite values of J\f. 

Figure 5 shows J\f (5) for the phenomenological model 
(3). The red line is the result of an average over the dis¬ 
tribution of the atoms in the mode while the inset with the 
blue line assumes maximally coupled atoms with fixed po¬ 
sitions. The averaged line shows a linear dependence with 
the Vacuum Rabi frequency that is proportional to y/N e s. 
We observe that the non-Markovianity increases with the 
number of atoms, and that there is a threshold for start¬ 
ing to observe non-Markovian behavior. This threshold can 
be understood intuitively: To have non-Markovian behav¬ 
ior the environment has to be correlated; to be correlated 
the environment has to consist at least of one atom. This 
threshold also bears an important message on the nature 
of non-Markovianity for open quantum systems: The pres¬ 
ence of oscillations in the system (non-zero vacuum Rabi 
oscillations) is not a good criterion for non-Markovian be¬ 
havior. Even in the non-averaged plot (inset of Fig. 5) 
the threshold at N ~ 1 to observe non-Markovianity is 
clearly visible. On the contrary, Rabi oscillations are al¬ 
ready present for N <C 1. 

Figure 5 also confirms that the reported speed-up is in¬ 
deed a non-Markovian effect since the speed of evolution 
grows with the number of atoms, which grows with the 
non-Markovianity [37]. Note that Fig. 5 also indicates that 
the non-Markovianity survives the averaging over the ran¬ 
dom distribution of atoms. 



FIG. 5. (color online) Measure of non-Markovianity (red line) in 
the cavity QED system (5) averaged over a distribution of atoms 
in the cavity mode as a function of the Vacuum Rabi frequency. 
The inset (blue line) shows the same A f for maximally coupled 
atoms. 


Concluding remarks . The reported experiment first 
demonstrated the theoretically predicted environment as- 
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sisted speed-up. Secondly considering the cavity field 
as a system coupled to a tunable environment realized a 
fully accessible quantum system obeying non-Markovian 
dynamics. Third demonstrated a novel way of environment 
engineering for optimal quantum control protocols [17]. 

The spatial distribution of the atoms, producing a broad¬ 
band of couplings, is playing a major role in allowing us to 
treat them as a non-Markovian reservoir as we do see the 
predicted behavior. 

These findings could prove useful to realize strong feed¬ 
back applied, e.g., in capture and release of a conditional 
state of cavity QED [38] and for a spontaneously created 
coherence [39]. In future work we will explore the quan¬ 
tum speed in optical cavity QED for large N to investigate 
any possible functional changes on the speed as suggested 
by Taddei et al. [16]. A measurement of the cross correla¬ 
tion between the field and the atomic polarization can give 
a complementary picture, and we are developing ways to 
do it. We finally note that the effect reported here is generic 
in the sense that it is not a peculiarity of two-level systems 
and we expect to go beyond it in forthcoming experiments. 
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